TU-745 
|hep-ph/0505252l 
May, 2005 



in ; Reconstructing Dark Matter Density with e + e Linear 

Collider in Focus-Point Supersymmetry 



(N 



(N '■ Takeo Moroi, Yasuhiro Shimizu and Akira Yotsuyanagi 

^ ■ Department of Physics, Tohoku University, Sendai 980-8578, JAPAN 

in 

(N 

in 
o 
in 
o 



Oh 



Abstract 



' It has been known that, in the focus point scenario of supersymmetry, the thermal 

relic of the lightest superparticle (LSP) is known to be a good candidate of the cold 



Q-i, dark matter. Assuming that the LSP in the focus-point scenario be the cold dark 

matter, we address a question how and how well the relic density of the LSP can be 



determined once the superparticles are found at future e + e~ linear collider. We will 
see that the determinations of the mass of the LSP as well as those of the Higgsino-like 
chargino and neutralinos, which will be possible by a study of the decay kinematics of 
the chargino or by threshold scan, will give us important information to theoretically 
reconstruct the relic density. Even if the Higgsino-like superparticles and the LSP are 
the only superparticles which are kinematically accessible, relic density of the LSP 
may be calculated with the accuracy of factor ~ 2; by adopting a mild theoretical 
assumption or by determining the masses of the Wino-like superparticles, uncertaintiy 
can be reduced to ~ 10 % or smaller. 



High energy physics and cosmology have had very important connections for deeper 
understandings in each field. In particular, progresses in physics at the energy frontier 
sometimes greatly improved our knowledges on the history of the universe. One of the 
most famous examples is the big-bang nucleosynthesis pQ; with the understandings of the 
interactions among the light elements and the expansion of the universe by general relativity, 
it became clear that the light nuclei were synthesized in the early universe. Remarkably, 
nowadays, light-element abundances can be precisely calculated and, for D, 4 He, and 7 Li, 
theoretical predictions are known to be in reasonable agreements with the observations [2]. 
As a result, we had an understanding of (some of) the components in the universe. This 
not only showed a very close connection between high energy physics and cosmology, but 
also provided a quantitative test of the big-bang scenario up to the cosmic temperature of 
~ 1 MeV. 

Now, we are in a position to have an understanding of physics at the electroweak scale. 
Although most of the results from the high energy experiments are more or less consistent 
with the predictions of the standard model of the particle physics, many particle physicists 
are expecting some new physics at the electroweak scale because of several problems in the 
standard model, like the naturalness problem, the hierarchy problem, and so on. Among var- 
ious possibilities, low-energy supersymmetry (SUSY) is the prominent candidate of physics 
beyond the standard model and hence signals from the superparticles are important target 
of the future high energy experiments. 

If the low-energy SUSY is realized in nature, it will also play important roles in cosmology. 
In particular, the lightest superparticle (LSP) is known to be a good candidate of the cold 
dark matter which, although there is no viable candidate within the standard-model particles, 
is strongly suggested by the present precise observations of the universe. Indeed, recent 
results from the WMAP suggests the dark matter density parameter to be [3J Ej 

(WMAP), 2 _ a ii O+0.008 /i N 

"DM n — u - iio -0.009' \ L ) 

where h here is the Hubble constant in units of 100 km/sec/Mpc. In the framework of the 
low-energy supersymmetry, however, it is not automatic to make the LSP to be the dark 
matter. In addition, although thermal relic of the LSP is the usual candidate of the dark 
matter, it has been also pointed out that the LSP dark matter may have non-thermal origin 
0. Thus, once the superparticles are found, it will be important to check if the thermal relic 
of the LSP is really the dominant component of the cold dark matter; if the theoretically 
calculated value of the LSP density well agrees with the observed dark-matter density, then 
it will be an important quantitative confirmation of the (simple) big-bang scenario up to the 
freeze-out temperature of the LSP, O(10 GeV). 

In order for a precise calculation of the relic density of the LSP, it is necessary to determine 
the properties of the superparticles by collider experiments. In particular, the relic density is 
sensitive to the properties of the LSP (such as the mass and couplings), precise determination 
of those parameters will be necessary. It has been discussed that the high energy e + e~ linear 
collider, recently named as the International Linear Collider (ILC), will be a good facility 
for a detailed study of the superparticles; for various cases, it has been shown that the e + e~ 
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linear collider can determine various mass and coupling parameters in the SUSY models 
very well (see, for example, Refs. [HI El IE])- Although, in the current situation, the LHC 
will probably be the first place where the superparticles will be found, it may not be easy 
to determine the SUSY parameters accurate enough to precisely calculate the relic density 
of the LSP P]. Thus, in the following, we consider how the ILC can help to determine the 
relic density of the LSP. 

Relic density of the LSP depends on how the LSPs annihilate in the early universe when 
they freeze out from the thermal bath. The dominant annihilation processes of the LSP are 
model-dependent and there are several parameter regions where the relic density of the LSP 
becomes consistent with the WMAP value, such as "bulk region," "coannihilation region," 
"rapid-annihilation funnels," "focus-point region," and so on ^U|. Thus, it is necessary to 
study the individual cases. 

In this letter, we consider one of the cases, the focus-point case; we assume that the 
dark matter consists of the LSP from the focus-point supersymmetric model, and address a 
question how and how well we can determine the density parameter of the LSP £\sp using 
the data from the ILC. With reasonable assumptions, we will see that, in the focus point 
case, relic density of the LSP can be well constrained. Before closing the introduction, we 
comment on how we proceed our analysis. We assume that the underlying model behind 
the minimal supersymmetric standard model (MSSM) is the focus point model (with grand 
unification) and that the dark matter is the thermal relic of the LSP. However, we will 
propose a procedure to calculate ^lsp without relying on any high-energy models (such as 
grand unification, mSUGRA-type parameterization, and so on); our aim is to determine all 
the relevant MSSM parameters from the collider experiments for the calculation of £\sp- 

We start with a brief review of the focus point scenario. The focus point scenario is 
characterized by the large (universal) scalar mass at the grand unification theory (GUT) scale 
^T].* 1 Although all the squarks and sleptons (as well as the heavier Higgs bosons) acquire 
multi-TeV masses, naturalness of the electroweak symmetry breaking may be maintained as 
far as the gauginos and Higgsinos remain relatively light. Thus, in this scenario, gauginos 
and Higgsinos are the "light" superparticles which may be accessible with the ILC with 
center-of-mass energy lower than ~ 1 TeV. 

Since the important superparticles are the charginos and neutralinos, the relevant param- 
eters for our study are the U(l)y and SU{2) L gaugino masses m G1 and mc2, the supersym- 
metric Higgs mass and the ratio of the vacuum expectation values of two Higgs bosons 
tan (3. Properties of the charginos and the neutralinos are determined by these parameters. 
In particular, the mass matrices of the charginos and the neutralinos are given by* 2 

M± = I -m G2 \ 
\ -y/2g 2 vsmf} \i H J 



^Such a large universal scalar mass may be a result of the "large cut-off supergravity" |12|. 

# 2 In fact, the chargino and neutralino masses may acquire radiative corrections, which should be taken into 
account in the actual study. In the focus point scenario, dominant radiative corrections are from the gauge 
boson loops and hence are calculable. In our study, we neglect the radiative corrections to the chargino and 
the neutralino masses. 
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/ —ma — g\v cos j3 

M _ -m G2 g 2 v cos (3 

— gxv cos (3 g 2 v cos (3 

V gtvsmp -g 2 vsm(3 fi H 

respectively. Here, g± and g 2 are the gauge coupling constants for U(l)y and SU(2) L , 
respectively, and v ~ 174 GeV is the (total) Higgs vacuum expectation value. These mass 
matrices are diagonalized by the unitary matrices U x + , U x - , and U x o as 

U X +M±U X - = diag(m x ±, "\±), U£ M.oU x o — diag(m x o, m x g, m x o, m x o). (4) 

If the GUT relation holds among the gaugino masses, the relation — 2mci holds (at 
the electroweak scale) while, for a successful electroweak symmetry breaking, becomes 
larger than |mci|. Then, the lightest neutralino, which becomes the LSP, is dominantly Bino 
although there is sizable contamination of the Higgsino, as we mention below. 

It should be noted that ^lsp is a lso determined by the previously mentioned four param- 
eters: m G i, m G2 , n Hl and tan/?. In the focus point scenario, the LSPs dominantly annihilate 
into the ti pair and into the gauge boson pairs {XiXi ~~ * W + W~ , Z°Z°) when the LSPs freeze 
out from the thermal bath ^HEl El- These processes are through the Higgsino components 
in x?; in the focus point scenario, can be relatively close to m G i and, consequently, 
sizable contamination of the Higgsino is possible in Xi- I* 1 this case, may become 

consistent with the WMAP value. In the cosmologically interesting parameter region where 
f^LSP ~ 0.1, \/j,h\ becomes smaller than |mc2| an d the lighter chargino Xi becomes Higgsino- 
like while the heavier one xf becomes Wino-like. In addition, in the neutralino sector, Xi 
(i.e., the LSP), x% X% an d xl are Bino-like, Higgsino- like, Higgsino-like, and Wino-like, 
respectively. As a result, the masses of xf > X21 an d X3 are quite degenerate with the approx- 
imate relations m x o ~ |^gi|? m x ± ~ m x ° 2 ~ m x° ~ \^h\i an d m x ± ~ m x l ~ l m G2|- These 
facts will become very important in the following study. 

As we mentioned, Qlsp is determined once the four parameters toqi, m G2 , fin, and tan/3 
are fixed. With the GUT relation, a narrow strip is obtained on the m G2 vs. fin plane, 
where ^lsp satisfies the WMAP constraint Eq.Q. Such a strip is insensitive to tan/5 (and 
the approximate relation [i H ~ 0.6m G2 holds on such a strip when m G2 ^ 300 GeV) [T4] . 
Although our strategy to determine ^lsp works for most of the points on the strip, (and 
even for cases without the GUT relation), it will be instructive to see the detail of several 
cases. Thus, we pick up two parameter points where WMAP value of fl-^sp is realized: Point 
1 with relatively small m G2 and Point 2 with larger m G2 . These points are listed in Table 
[U for these points, f2 L sp an d the lightest Higgs mass m h are calculated with adopting that 
all the sfermion masses are 3 TeV. In the following, several physical quantities such as the 
cross sections as well as the estimated errors in the reconstructed f^gp ~vvill be given for these 
points. (In our numerical analysis, we use the DarkSUSY package [T3] to calculate ^lsp-) 

Now we discuss the role of the ILC. As we mentioned, in the focus point case, the 
chargino(s) and neutralino(s) are the ones which can be produced and studied at the ILC. 
Thus, we focus on the question what can be measured and studied by the production of the 
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Point 1 Point 2 



togi 144 GeV 240 GeV 

m G2 300 GeV 500 GeV 

Hh 200 GeV 307 GeV 

tan/3 10 10 

m x o 127 GeV 226 GeV 

m x o 190 GeV 304 GeV 

m x o 208 GeV 311 GeV 

m x J 335 GeV 522 GeV 

m x ± 176 GeV 291 GeV 

ml 335 GeV 521 GeV 

m h 116 GeV 116 GeV 

n LSP h 2 0.113 0.113 

Table 1: Points to be used for our analysis. 



charginos and the neutralinos.* 3 For this purpose, we first calculate the production cross 
sections of the charginos and neutralinos. Since the sleptons are extremely heavy, chargino 
and neutralino production processes are mediated by the s-channel gauge-boson exchange 
diagrams; neglecting the selectron diagrams, the chargino and neutralino production cross 
sections are given by 



XxXy) = 



8tt 



Cll\ 2 + \Clr\ 2 



E X E Y + - |p f | 2 



+ (CllC* lr + C* LL C LR ) m x m Y 



(5) 



where Ex, Ey, and pf are the energies of xx and xy, and their three- momentum, respec- 
tively, N X y = \ when xx and xy are identical and N X y = 1 otherwise, and 



s 2 — 2 (m x + m Y ^j s + (m x — m Y ^j 



(6) 



with yfs being the center-of-mass energy of the ILC. Here, the polarization of the electron is 
specified to be left-handed while the average over the positron polarization is taken. Result 
for the right-handed electron is given by a similar formula by replacing Cll —> Crl and 
Clr — ► Crr. For the chargino production, 



[Cll\ y ± y =f 



— $XY + 




12X 



X \2Y 



(7) 



#3 We assume that the squarks and sleptons will be known to be very heavy by the study of the LHC, and 
also by the negative searches of these particles at the ILC. 
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Figure 1: Production cross sections of the charginos and neutralinos for Point 1 for the 
processes as functions of y/s. Here, we have neglected the selectron-exchange diagrams and 
averaged over the polarization of the electron beam. 



where g\ = g\ + gr|, and e = (g 1 2 + g 2 2 ) 1//2 is the electric charge. [Crl\ x ± x t is obtained by 
replacing g e z L — > g e z R , where 



n e L _ 9l flf e R _ 9l 

9z - -TT > 9z - —■ 
z 9z 9z 



In addition, \Clr) v ± „=r and \Crr\ + i are obtained by the following replacements: 

X X X Y J AxXy 



\- Clr ^x x x y 



^ Cll ^x x x y 



\- Cr ^x x x y 



\- Crl ^x%x y 



For the neutralino production processes, 



[Cll] x xX y 



9z L 9z 



2(s 



VTir 



l U X°]lx \ U X°\zY ~ \ U X°\lx \ U X°- 



AY 



(9) 



(10) 



and [C lr } x o xX o y 
replacing g z L — 



- [C L l] XyXx - 



For the right-polarized electron, the results are obtained by 



9z 



To see the size of the cross sections, we plot the cross sections for various final states 
in Figs. ^ and |2 for Point 1 and Point 2, respectively, as functions of yfs. As one can see, 
for the processes e + e~ — > XiXi an d X2X3; cross sections become significantly large (if these 
processes are kinematically allowed). As we will see, these processes give very important 
information in calculating the relic density of the LSP.* 4 On the contrary, cross sections 

# 4 Notice that the chargino \t mostly decays into W^^Xi final state (where the superscript * is for 
virtual particles) while %° an d X3 decay into Z°^*'x\ an d ^X? final states (where h is the lightest Higgs 
boson). Thus, we expect that some of the chargino and neutralino production events can be distinguished 
by selecting particular final states. For example, for e + e~ — > X1X1 > we can use 11 + 2 jets + missing E^ 
events, which is not from X2X3 production, and for e + e~ — > X2X3J 2 I + 2 jets + missing Et events can be 
used. 
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Figure 2: Same as Fig. ^ except for Point 2. 

cr(e + e~ — ► X2X2) an d &(e + e~ — > X3X3) are ver y small. This is due to the facts that the cross 
sections acquire suppressions from the mixing factors and that these are p-wave processes. 
Note also that the measurements of the cross sections for e + e~ — ► XiXi an d X2X3 wu l tell 
us that Xi ) X2i an d X3 are Higgsino-like, not Wino-like. 

A great advantage of the ILC is that the beam energy can be easily tuned. Thus, it is 
possible to study various production processes at the threshold region. In particular, in the 
case of the focus-point scenario, threshold scans of the chargino and neutralino productions 
give us important information. From the chargino production process e + e~ — *> X1X11 we 
can have a precise measurement of the (lighter) chargino mass Tn x ±. Study at the region 
y/s ~ m x o + m x o is also interesting. At y/s ~ 2|/i#|, three different neutralino production 
processes become kinematically allowed since the masses of X2 an d X3 are close. However, 
the cross sections for e + e~ — > X2X2 an( ^ e + e~ — > X3X3 are very small in particular at the 
threshold region since these processes are p-wave suppressed; we can see that the cross 
sections for these processes are much smaller than a(e + e~ — > xtXs)- Thus, we observe only 
the process e + e~ — > X2.X3- Since the process e + e _ — > X2X3 becomes kinematically allowed 
when yfs = m x o + m x o, we can determine the combination m x o + m x o or, equivalently, the 
averaged mass of X2 an( i 

m x° 23 = 2 { m x°2 + m xl) > ( U ) 

by the scan around y/s ~ m x o + m x o . 

In order to study the properties of the LSP, on the contrary, y/s should better be op- 
timized so that the production of the lighter chargino, which decays into the LSP Xii i s 
enhanced. Once the charginos are copiously produced, then the mass of the LSP (more pre- 
cisely, the mass difference m ± — m x o) is determined by the study of the energy distribution 
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of the decay products.* 5 

At the ILC, errors in the measurements of m x ±, ^ x o 3 , and m x o are expected to be mostly 
from the detector resolutions [E]. For example, it was pointed out that, for some choice of 
the SUSY parameters, masses of the charginos can be determined using e + e~ colliders with 
the errors of ~ 50 MeV by the threshold scan. In addition, from the energy distribution of 
the decay products of the chargino and neutralinos, the mass of the LSP is also determined 
with the uncertainty of ~ 50 MeV. Although these results are for the case of Wino-like 
chargino and neutralino, we expect that three mass parameters (i.e., Tn x ±, m x o, and ^ x o 3 ) 
are accurately measured once xf-> X21 an< ^ X3 become kinematically accessible at the ILC. 
Since Xi is Bino-like while xt ( as we U as X2 anc ^ X3) are Higgsino-like, we can constrain moi 
and hh from the measurements of m x o and m x ± (or from the masses of other Higgsino-like 
neutralinos). In addition, the "mass difference" rh x o 3 ~ m x ± is sensitive to some combination 
of tan (3 and mc2- 

Although the relic density of the LSP depends on four MSSM parameters (mci, mc2, 
/Jijj, and tan/3), interesting bound on Qlsp can be obtained even at this stage. To see this, 
we can perform the following analysis. Let us imagine a situation where fn x ±, m x °j anc ^ 
m x o a are well measured at the ILC. Using these quantities, we impose three constraints on 
the four underlying parameters and determine tugi, ^h, a nd tan/5 as functions of m G 2- In 
the determination of mqi and in fact, there are four possible choices of their signs: 
(sign(m G i), sign(/i#)) = (+, +), (+, — ), (— , +), and (-, -). #6 Effects of the signs of fi H and 
tugi are quite different. In order to see how the reconstructed relic density depends on m G2 
and sign(//#), here let us consider the case where the sign of the reconstructed m G1 is the 
same as that of the underlying one; effects of sign(m G i) will be discussed later. Once we 
reconstruct mci, /jlh, an d tan/3, we calculate the relic density of the LSP as a function of 
iriQ2i which we call 

^lsp (m G2 ; m x ± , m x o , rh^ ) . 

In Figs. El and El we plot ^LSp( m G2; m x ±? m x?' ™xt?) as a f unc tion of m G2 , with Tn x ±, 
m x o, and ^ x o 3 being fixed to be the values from Point 1 and Point 2, respectively. We use a 

fixed value of the lightest Higgs boson mass m^. (We have checked that ^lsp is insensitive to 
the change of rrih-) The lines have endpoints; this is due to the fact that, when m G 2 becomes 
too large or too small, there is no value of tan (3 which consistently reproduces the observed 
mass spectrum. To demonstrate this, we also showed the points where tan/3 takes several 
specific values. Here, one should note that there exists two-fold ambiguity from the sign of 
fin] in our study, we reconstructed \Xh for two cases (with the sign of the underlying value 

#5p rom t^-O^o production, we may in principle perform a similar analysis to determine the mass of the 
LSP. In this case, however, error may become larger because the masses of X2 an d X3 differ by 0(1 GeV). 
Global fit using all the data should be appropriate once the charginos and neutralinos are really found in 
the future. 

^ 6 To be more precise, these signs are the relative signs between mci and mG2 or /ih and rriG2- We assume 
that the gaugino masses and fin are real in order to avoid constraints from CP violations. 
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Figure 3: ^lsp { m G2 ] m x ± ? m x\ > ™x° 3 ) as a f unc ti° n of m G2) where tti x ±, m x o, and ^ x o g are 
fixed by the underlying values for Point 1 with positive fxn (solid) and negative /xh (dashed). 
Marks on the figure indicate the points with tan f3 — 2, 3, 5, 10, 50. 



of ix h being fixed). In the figures, the results for the cases with positive and negative hh are 
shown. As one can see, the line for the fin < case is "attached" to one of the endpoints of 
the line for llh > 0. This is from the fact that, at the tree level, flip of sign(/i#) is equivalent 
to the change (5 — > n — (5. 

As one can see, the case with negative hh may give large uncertainty in the reconstructed 
f2 LSP . If ix h < 0, however, smaller wq2 is required than in the case of positive-/x# in order 
to reproduce the observed mass spectrum, as shown in the figures. Then, m x ± , for example, 
may become smaller than the experimental bound from the negative search for the X1X2 
production process. For Point 1 (Point 2), ^^480 GeV (750 GeV) is enough to exclude 
the hh < case. Thus, in the following, we assume that this is the case and neglect the 
uncertainty from the Hh < case. Then, even without any further constraint on m G2 , relic 
density of the LSP can be determined within a factor of ~ 2 or smaller. 

For a quantitative study of the uncertainties in the reconstructed LSP density, we define 
the following quantity: 



d In fi LS p 



dlnmi 



mi = m xf ' m x\ ' m x§ 3 ' m G2 , 



(12) 



so that* 7 
<^lsp _ 

^LSP 



D 



m ± 



^ m xt 



+ D 



6m x\ 

m n 

1 m x° , 



+ D 



X23 



m n - 
x 23 17l v 

X23 



+ D 



5m, 



G2 



m G2 



m G2 



1/2 



.(13) 



# 7 In fact, some of the observables may be correlated (in particular because, as we mentioned, one of the 
observable would be m ± — m o, not to o). Here, we neglect effects of such correlations. 
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Figure 4: Same as Fig. El except for Point 2. 
For Point 1 and Point 2, we found (D m ± ,D m ,D m ,D mG2 ) to be (19,8.1,8.6,1.5) and 

x l x l x 23 

(7.5, 7.9, 8.5, ~ 0.01), respectively. (Here, notice that D mQ2 for Point 2 is accidentally small 
because of the choice of the underlying parameters.) As we mentioned, at the ILC, m ±, 
m x o, and m x o 3 will be well determined with the error of 0.1 % or less jHj (although the errors 
should depend on the model parameters). Using the values of D m ± , D m , and given 

x l x l x 23 

above, <K)lsp from 5m x ±, 5m x o, and St?i x o 3 is very small. 

For a better determination of I^lsp, some information about mo2 is necessary. For exam- 
ple, if we try to determine £)lsp with ~ 10 % accuracy, which is the level of the WMAP value, 
mQ2 should be determined within the uncertainty of ~ 30 GeV for Point 1 and ~ 170 GeV 
for Point 2. If we impose the GUT relation among the gaugino masses, mQ 2 is accurately 
determined. However, we do not pursue this direction since we hope to calculate Qlsp in 
the framework of the MSSM. 

The SU(2)l gaugino mass mc2 is approximately equal to the masses of xf an d X% Thus, 
if we know something about m ± or about m Y o, it will help us to set a bound on m G2 . Some 
constraint may be given by the LHC. Even in the focus-point scenario, a sizable amount of 
the gluino may be produced at the LHC, and some of them decay into xf or X°i (and other 
quark jets). These charginos and neutralinos will decay by emitting W ± or Z° boson (or the 
lightest Higgs boson, in some case). If the gauge bosons subsequently decay into the leptons, 
for example, we may have events with multi-leptons and missing whose background is 
probably relatively small. By studying the energy distribution of the charged leptons, some 
information about the Wino-like particles may be obtained. 

If xf or Xa becomes kinematically accessible at the ILC, some direct information of these 
particles will become available. For example, we can look for the associated production 
processes of these Wino-like particles with lighter charginos or neutralinos (although the 
cross sections for the associated productions are relatively suppressed). For example, with 
large enough y/s, cross section for the xt X2 production process is sizable, as shown in Figs. 




9 



H] and El With the threshold scan of this process, m x ± may be measured, resulting in a good 
determination of mQ 2 . 

In order to constrain mc2, we may also use tan/3 dependence of the lightest Higgs mass. 
As shown in Figs. El and EJ the reconstructed value of tan/3 depends on mc2- Importantly, 
smaller value of tan (3 results in more suppressed value of the lightest Higgs mass. Although 
the lightest Higgs mass depends on other MSSM parameters (in particular, on the stop 
masses [IE]), too small tan/3 will be excluded by the detailed study of the lightest Higgs, 
which will be possible at the ILC [HI 13 IE]- If we obtain a lower bound on tan/3 from the 
measurement of the Higgs mass, it will also help to determine the relic density of the LSP. 

Another possibility to constrain mc2 is to measure the cross sections for the processes 
e + e~ — > X1X2 an d e + e~ — > XiXz- ^ n some case, X2 an d X3 both dominantly decay into Z°(*'Xi- 
If so, the X1X2 an d X1X3 production processes has the final state with 2 / + missing and 
2 jets + missing E T . It may not be easy to distinguish these two events, but we can just 
consider the total cross section <r(e + e~ — > X1X2) + cr(e + e~ — > xIxT)- F° r Point 2, for 
example, the total cross section a(e + e~ — > X1X2) +o r (e + e~ — > X1X3) monotonically increases 
as a function of 777,02 from 15 fb (for mc2 = 440 GeV) to 29 fb (for 777-02 = 1830 GeV). 
Thus, if this cross section is precisely measured, we can obtain another information about 
mG2- However, it should be noted that several serious standard-model backgrounds may 
exist. In particular, the processes e + e~ — > W + W~, u e u e Z°, and W + W~Z° have large cross 
sections. (For ^1 = 1 TeV, a(e + e~ -> W + W~) ~ 3 pb, a(e + e~ -> v e v e Z Q ) ~ 900 fb, 
and cr(e + e _ — > W + W~ Z°) ~ 60 fb for the unpolarized electron beam |T7j.) Although 
these processes are suppressed for the right-polarized electron beam, cross sections for these 
processes are much larger than er(e + e~ — > X1X2) +0"(e + e~ — > XiXf)- Thus, some appropriate 
kinematical cuts should be necessary to use these modes for the determination of <r(e + e~ — > 

Xlxl) + <y{e + e~ ^ xlxl)- 

So far, we have not considered effects of the sign of m G1 . Unfortunately, f2 L sp depends 
on the sign of mci although the determination of sign(mci) seems challenging. For the case 
with negative mci, ^lsp^ 2 varies from 0.9 to 7.4 (Point 1) and from 0.2 to 6.7 (Point 2). 
If sign (mci) is undetermined, thus, two- fold ambiguity will remain. However, experimental 
determination of sign(mQi) may be possible [TS|. In addition, if the GUT relation among 
the (absolute values of) gaugino masses is experimentally confirmed, it will give another hint 
of the signs of the gaugino masses. 

We have not discussed possible errors of £\sp originating from parameters other than 
m x ±, m x o, and ^ x o 3 , and mc2- Theoretical prediction on ^lsp, in fact, depends also on 
other parameters. For the process X1X1 ~^ ^i, which can be the dominant pair annihilation 
process of the LSP, there exists the t channel stop exchange diagram, so the cross section for 
this process depends on the stop masses. Since the stops are very heavy in the focus point 
scenario, the stop masses will not be directly measured at the ILC. In addition, the cross 
section for this process also depends on the lightest Higgs boson mass through the s-channel 
exchange of the Higgs boson. (Of course, the lightest Higgs boson mass will be accurately 
measured at the ILC, so the latter diagram will not give serious uncertainty.) However, the 
process X1X1 ~^ tt is dominated by the s-channel Z° exchange diagram 0. We have checked 
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that the reconstructed value of £\sp is insensitive to the choice of the stop mass and m^. 
Other possible uncertainty may be from the pair annihilation into the bb final state, which 
may give ~ 10 % contribution to the total annihilation cross section when tan/3 is very as 
large as ~ 50 because of the enhanced bottom-quark Yukawa coupling The cross section 
for the process X1X1 ~ > bb cannot be calculated unless we know the pseudo-scalar Higgs 
boson mass. Thus, ~ 10% uncertainity may remain unless the pseudo-scalar Higgs boson 
mass is somehow constrained. Uncertainties from the first and second generation sfermions 
are irrelevant since the pair annihilation into light fermions is p-wave suppressed. 

In this letter, we have seen that, if the dominant component of the cold dark matter 
is the thermal relics of the LSP in the focus-point supersymmetric model, we will have a 
good chance to theoretically reconstruct the cold dark matter density once the superparticles 
becomes kinematically accessible at the ILC Positive confirmation of the dark matter density 
will give us a quantitative understanding of our universe up to the freeze-out temperature 
of the LSP, O(10 GeV). Of course, even in the framework of the supersymmetric models, 
there are several possibilities to realize the LSP dark matter. Superparticle spectrum and 
the observables at the collider experiments depend on the model. Thus, for other cases, it 
is also necessary to study how and how well we can reconstruct the dark matter density. In 
addition, in our study, we have not performed detailed detector simulations to have accurate 
errors in the observables (although we have seen that, using the reasonable estimate of the 
errors, determination of ^lsp seems promising with a good accuracy). We leave these for 
future studies. 

Note added: While finalizing this letter, we found the paper J§| which may have some 
relevance with our analysis. 
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